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Abstract 

We discuss jets in strongly coupled N — 4 supersymmetric Yang-Mills plasma and their dual 
gravitational description. 



1. Introduction 

The discovery that the quark-gluon plasma produced at RHIC behaves as a nearly ideal fluid 
|1 1 has prompted much interest into the dynamics of strongly coupled plasmas. Energetic par- 
tons produced in hard processes during the early stages of heavy ion collisions can traverse the 
resulting fireball and deposit their energy and momentum into the medium. Analysis of particle 
correlations in produced jets can provide useful information about the dynamics of the plasma 
including the rates of energy loss and momentum broadening J2]|3], as well as the speed and 
attenuation length of sound waves |3J. 

Studying the dynamics of jets in QCD from first principles — from the initial hard process 
which created the jet to energy loss, showering and thermalization — is a challenging task. Be- 
cause of the lack of theoretical tools available for analyzing strongly coupled dynamics in QCD, 
it is useful to have a toy model which is qualitatively similar to QCD but where strongly coupled 
dynamics can be treated in a controlled manner. A class of such toy models are non-Abelian 
theories with gravitational duals 0[5]. The most widely studied example is that of strongly cou- 
pled N = 4 supersymmetric Yang-Mills theory (SYM). The deconfined plasma phases of QCD 
and SYM share many qualitative properties. For example, both theories describe non-Abelian 
plasmas with Debye screening, finite spatial correlation lengths, and long distance dynamics de- 
scribed by neutral fluid hydrodynamics. When both theories are weakly coupled, appropriate 
comparisons of a variety of observables show rather good agreement (6j [71 [8] . This success, 
combined with the lack of alternative techniques for studying dynamical properties of QCD at 
temperatures where the plasma is strongly coupled, has motivated much interest in using strongly 
coupled SYM plasma as a model for QCD plasma at temperatures T of a few times Aqcd (or 
1.5 T c < T <4T C ). 

Despite their qualitative similarities at finite temperature, it should be emphasized that dy- 
namics in SYM are not a controlled approximation to those in QCD. Moreover, in regimes where 
asymptotic freedom may be important the qualitative dynamics in SYM may not resemble those 
in QCD. This is certainly the case during the initial hard process which creates a jet. In QCD, 
asymptotic freedom guarantees that for sufficiently high energy processes, the early time dy- 
namics will be weakly coupled. Therefore, in QCD the production of high energy jets will be 
governed by perturbative physics. This should be contrasted with SYM, where the coupling 
doesn't run and early time dynamics are just as strongly coupled as late time dynamics. 
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Because of the above point, we chose to focus of quantities which are largely insensitive to 
the initial event which created the jet. Our motivation is similar to that of Ref. (9), in which 
weak coupling physics in asymptotically free QCD is envisioned as producing a high energy 
excitation, whose propagation through the plasma is then modeled by studying the behavior of 
a similar type of excitation in a strongly coupled SYM plasma. As SYM contains only adjoint 
matter, fundamental fields must be added to the theory in order to study the propagation of 
fundamental representation excitations. Within the context of gauge/gravity duality, a natural 
choice is to add a fundamental N — 2 hypermultiplet to SYM and study N-2 excitations 1 10]. 

By studying N-2 excitations in strongly coupled SYM plasma we can then answer the fol- 
lowing simple set of important questions. How far can an energetic jet in a plasma propagate 
before it thermalizes? Where along its trajectory does a jet prefer to deposit energy and mo- 
mentum? How well is the radiation created by the jet modeled by hydrodynamics? Answering 
this last question is important as hydrodynamics is a much simpler theory to work with than any 
complete microscopic theory and it applies equally well to QCD and SYM. 



2. Jets and their dual gravitational description 

Energetic quarks moving through a plasma are quasi-particles — they have a finite lifetime 
which can be long compared to the inverse of their energy. Because of this, some care is needed 
in defining observables such as the instantaneous energy loss rate or the penetration depth. At 
weak coupling, where a perturbative analysis is applicable, an energetic quark scattering off 
excitations in the medium can emit gluons which may subsequently split into further gluons or 
quark-antiquark pairs. An energetic quark may also annihilate with an antiquark in the medium. 
A natural question to consider then is which quark should one follow when computing either the 
position of the total excitation or its instantaneous energy loss rate? Once a quark has emitted 
a qq pair, or annihilated with an antiquark, it becomes ambiguous which quark was the original 
one. This issue is most cleanly avoided if one focuses attention not on some (ill-defined) "bare 
quark", but rather on the locally conserved charges of the entire dressed excitation. 

In both QCD and SYM coupled to a fundamental N-2 hypermultiplet, there are two con- 
served quantities which are useful for characterizing a fundamental representation excitation. 
These are simply the baryon number (or quark number) current a and the total stress tensor 
for the system T^ v . Even though qq pairs can be produced by an energetic quark traversing the 
plasma, conservation of /£ aryon and T MV implies that the total baryon number and four momentum 
of the state will remain constant. For an energetic, well collimated jet, the baryon density can 
remain highly localized for a long period of time. It is the collective excitation with localized 
baryon density which we will refer to as a dressed quark, or for the sake of brevity, simply as a 
quark. 

With the above points in mind, let us for the moment focus not on the microscopic details 
needed to fully describe a quark moving through a plasma but rather on the large scale structure. 
Fig. [T] shows a cartoon depicting the large scale picture of a quark moving though a plasma. 
The quark itself is a localized object whose location is specified by the localized baryon density. 
Specifically, the trajectory of the quark is defined to be 

*quark(0 = ^ J d 3 X X p(t, X), (1) 

where p = J®. dryon and Q = f d 3 xp(t,x) is the baryon number of the quark. For very energetic 
quarks the stress tensor will also be highly peaked in the vicinity of JC quar k(0- As the quark 
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Figure 1 : A cartoon depicting the large scale structure of an energetic quark moving through a plasma. 



moves through the plasma, it will experience frictional drag forces and consequently, energy 
and momentum will be transported from the vicinity of jc qua rk(0 to the far zone. Because SYM 
is dissipative, gradients in the stress tensor necessarily become smaller and smaller at larger 
distances, and the far-zone dynamics of emitted radiation are governed by hydrodynamics. In 
the long wavelength limit the total stress tensor may be written 



JVUV 



hydro quark' 



(2) 



where T^ ik is localized about x qua rk(0 and T^ dm satisfies the constituent relations of hydrody- 
namics together with the energy-momentum conservation equation 



d u T\ 



hydro/*' f ) = f )' JV ( X < f ) S -^Crk^' 



(3) 



The effective source P for hydrodynamics is localized at the position of the quark and takes 
into account all microscopic dynamics in the near-zone relevant for the determination of the far- 
zone hydrodynamic behavior. In general its determination will require a complete microscopic 
calculation. However, just by energy-momentum conservation, one crucial property the effective 
source must satisfy is that its integral over space must give the rate that the quark loses four- 
momentum: 



d 
dt 



(4) 



where fif) is minus the drag force acting on the quark, and thus the rate that four-momentum is 
transfered to hydrodynamic modesQOperationally, the instantaneous drag force can be obtained 
by enclosing the quark's baryon density in a sphere and computing the total outward flux of four 
momentum as a function of time. 

Just as in electromagnetism, where local sources can be expanded in terms of point-like 
multipole moments (i.e. derivatives of delta functions) with each moment producing a particular 
falloff of the fields with distance, the long distance hydrodynamic response of an energetic quark 
traversing a plasma can be systematically expanded in inverse powers of distance. In doing so, 



1 More precisely, as hydrodynamics is an effective theory valid on scales much larger than the mean free path, the 
integral of the effective source must only give the microscopically time-averaged drag force acting on the quark. 
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the leading long distance behavior of the hydrodynamic response is determined by the source's 
monopole moment IfTTI . More explicitly, at leading order the source may be approximated 

nx,t)^f(t)6\x-x quaik (t)). (5) 

We note that the coefficient of the delta function is constrained via Eq. Q. Together with the 
hydrodynamic constituent relations, Eqs. Q and <|3j show that the leading order behavior of the 
far-zone stress tensor is determined solely by the quark's trajectory Jr qua ,k, f and the plasma's 
viscosity. 

In general there are two possible hydrodynamic modes an energetic quark can excite. These 
modes are a sound mode and diffusion mode. The diffusion mode trails behind the quark and 
slowly broadens with distance. This mode, symbolized in Fig. [T]by the trailing flow lines, con- 
sists of laminar flow in the same direction as the quark. For supersonic motion the sound mode 
manifests itself as a Mach cone or sonic boom. This mode is symbolized in Fig.[T]by the outward 
moving conical wake emanating from the quark. 

For heavy quarks whose motion is not ultrarelativistic, f(t) only depends on the temperature 
of the plasma, the 't Hooft coupling and the quark's velocity lTT2l fT3l IT4TI . However for ultrarela- 
tivistic quarks f(t) can be sensitive to the initial conditions which created the quark fl5l . As the 
weakly coupled hard process which creates a jet in QCD is very different than that of strongly 
coupled S YM, for the case of light quarks — whose motion is always ultrarelativistic — it is de- 
sirable to have another observable which is related to energy loss but not very sensitive to initial 
conditions. One quantity to consider is the maximum distance Ax max (E) which a quark with ini- 
tial energy E can travel. One can imagine fixing the quark's energy and varying all other degrees 
of freedom characterizing the state and looking for the particular configuration which maximizes 
the penetration depth. Therefore, essentially by construction, the initial conditions associated 
with quarks that maximize the penetration depth at fixed energy are much more constrained. 

To define the penetration depth, we imagine measuring x qu!1 ±(t) at some early time f*. We 
then define the penetration depth Ax as Ax = |jr qU ark( 00 ) - x qud ±(t t )\ . To make the quantity Ax 
meaningful, we also need to measure the quark's energy at time ?*. If the quark has been moving 
for some time prior to t t , it will have deposited energy into the plasma — we must disentangle the 
energy deposited in the plasma from the remaining energy of the quark itself. In the limit where 
the quark has an arbitrarily large energy which is arbitrarily localized about x quai k, separating the 
quark's energy from the energy transfered to the plasma will be unambiguous. We consider only 
this limit in this paper. 

We now discuss jets in SYM from the 5d gravitational perspective. The lOd geometry cor- 
responding to an infinite equilibrium SYM plasma is the product of a 5d sphere and the 5d 
AdS-Schwarzschild (AdS-BH) geometry, whose metric is ds 1 = ^ [-f(u) dt 2 + dx 2 + du 2 /f(u)], 
where f(u) = l-(u/u/,) 4 and L is the AdS curvature radius. The AdS-BH geometry contains a Ad 
boundary at radial coordinate u — 0. This boundary has the geometry of Minkowski space and 
is where the Ad dual field theory should be thought of as residing. The AdS-BH geometry also 
contains a translationally invariant event horizon at u = m/, and has a corresponding Hawking 
temperature T = l/nu/,. 

The addition of a N — 2 hypermultiplet to SYM is accomplished in the gravitational setting 
by adding a D7 brane to the KM geometry IfTUl . The D7 brane fills a volume of the AdS-BH 
geometry which extends from the boundary at u — down to maximal radial coordinate u m . The 
bare mass M of the hypermultiplet is proportional to l/u m lfT2l . Therefore, for massless quarks 
the D7 brane fills all of the five-dimensional AdS-BH geometry and for very heavy quarks the 
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D7 brane ends very close to the boundary at u = 0. Open strings which end on the D7 brane 
correspond to dressed quarks in the dual field theory. 

In general, to study the dynamics of qq pairs in the field theory one must study dynamics 
of strings in the full lOd geometry. The reason for this is that in general string endpoints have 
non-trivial trajectories on the S 5 . However, in the limit of large quark mass and very small quark 
mass it is consistent to set the motion on the S 5 to vanish [15|. Because of this simplification 
and because any additional motion of the string on the S 5 will only add to the energy of the 
string without affecting its stopping distance in the spatial direction, we choose to study only the 
dynamics of very heavy and massless quarks in this paper. 

The motion of strings corresponding to very heavy or very light quarks is quite different. For 
the case of very heavy quarks, the D7 brane ends at u — > and string endpoints always remain 
close to the boundary. For very light quarks, the D7 brane fills the entire AdS-BH geometry 
and string endpoints are allowed to fall unimpeded towards the black hole 03J. However, in 
either case the late-time spatial motion of the string endpoint always ceases and in either case 
the total distance traveled by the endpoint can be made arbitrarily far by giving the string more 
momentum in the spatial directions IT5l . 

Fig. [2] shows a cartoon illustrating the gravitational dual of a quark moving through a SYM 
plasma. This consists of a string (shown as the magenta curve) moving in the AdS-BH geometry. 
The endpoints of strings are charged under a t/(l) gauge field whose influence extends to the 4d 
boundary. The boundary itself behaves as an ideal electromagnetic conductor [ 16 1, so the pres- 
ence of a string endpoint induces an image current density on the boundary above the endpoint. 
This is illustrated schematically in Fig. [2] Via the standard gauge/gravity dictionary [4, 5| the 
induced current density corresponding to the string endpoint has the field theory interpretation as 
minus the expectation value of the baryon current density of a dressed quark. In the limit where 
the string endpoint is close to the boundary the induced baryon density is highly localized and 
we may approximate x quai k ~ x str i„<,, where x stI ing is the spatial position of the string endpoint 
03. 




Horizon 



Figure 2: A cartoon illustrating the gravitational dual of a quark moving through a SYM plasma. 



Via Einstein's equations, the presence of the string will also perturb the 5d geometry. In the 
large N c limit, the corresponding perturbation is order 1 /N% and may be analyzed by linearizing 
Einstein's equations. As in the electromagnetic problem, classical field theory implies the near- 
boundary perturbation in the 5d geometry will induce a corresponding perturbation in the Ad 
boundary stress tensor lfl8l[T9ll . Via the gauge/gravity dictionary, the stress tensor induced on 
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the boundary has the interpretation as the expectation value of the complete field theoretic stress 
tensor, and hence contains information about dynamics on all length scales. 

If the string endpoint is close to the boundary — either because of initial conditions or the 
restriction of a heavy quark — the string endpoint will have a high gravitational potential energy. 
It is this UV sensitive part of the string energy that we identify with the energy of a quark. Via 
the gravitational equivalent of the electromagnetic image problem (i.e. the gravitational bulk- 
to-boundary problem), the high energy and momentum densities near the string endpoint get 
mapped onto a region of Ad space which roughly coincides with the location of the quark's baryon 
density. That is, 7^ ark , as defined above in (Eb, is induced in the gravitational framework by the 
high gravitational and kinetic energy carriea by the string endpoint. Again, this is illustrated 
schematically in Fig. [2] 

With the above point in mind, a natural question arises: what is the 5d gravitational descrip- 
tion of the drag force f**l The answer to this question can be understood by noting that each 
conserved "charge" of the string gets mapped into a conserved charge in the field theory. In the 
gravitational description, there in a conserved flux of energy and momentum which flows down 
the string from the endpoint. As the energy and momentum carried by the near-endpoint portion 
of the string induces T^^, it must be that that this flux is related to f. In the long wavelength 
limit, retardation effects in the propagation of gravitational disturbances from the bulk to the 
boundary can be neglected and the locally conserved flux of energy and momentum down the 
string — averaged over "microscopic" time scales of order the horizon radius — must in fact 
coincide with 

3. Results and discussion 

Heavy quarks decelerate slowly. Therefore, as long as one focuses on local questions, for 
sufficiently large mass it is useful to consider the limit where the quark's velocity v is constant. 
Via gauge/gravity duality, the energy loss rate for a heavy quark moving through SYM plasma at 
constant velocity evaluates to dE/dt = —it ^[AT 2 v 2 1(2 Vl - v 2 ), where A is the 't Hooft coupling 
ifTZl [T3l fl4ll . The complete (i.e. valid on all length scales) energy density and energy flux of a 
quark moving at constant velocity has been evaluated in Refs. |[TTll20ll2Tll22l . Fig. [3] taken from 
Ref. [11], shows the complete energy density and energy flux created by a heavy quark moving 
at velocity v = 3/4 though a strongly coupled SYM plasma. In a conformal theory such as SYM, 
the speed of sound is c s = VI /3 ~ 0.58, so Fig. pi shows supersonic motion. As is evident 
in Fig. [3] a Mach cone, with opening half angle 8m as 50° is clearly visible in both the energy 
density and the energy flux. Near the Mach cone, the bulk of the energy flux flow is orthogonal 
to the wavefront. Also evident in Fig. |3]is the presence of a diffusion wake in the energy flux. 

As emphasized in the previous section, with knowledge of the drag force on the quark and 
the viscosity of SYM plasma, one can completely determine the large-distance structure of the 
quark's wake. A quantitative comparison of the complete quark wake to that predicted by hydro- 
dynamics with the leading order effective source <(3j was done in Ref. fTD . and it was demon- 
strated that the hydrodynamic approximation provides a good description of the quark wake even 



2 This statement requires further explanation. Hydrodynamics is valid on scales much greater than the microscopic 
scale l/T. Retardation effects in the bulk to boundary problem are order i/T and hence to leading order can be neglected 
in the hydrodynamic limit. As noted in footnote[T] when one extracts f(t) from hydrodynamic sources, as in Eq. |4j, 
one only extracts the drag force averaged over microscopic times. Therefore, one should regard the energy-momentum 
flux down the string as being equivalent to ff(f) only up to course graining in time. 
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Figure 3: Left— Plot of energy density |jc|A£(jc)/(T 3 VI). Right— Plot of energy flux \x\AS(x)/(T 3 VI). In both plots 
the quark is moving in the x\\ direction at speed v = 3/4. The flow lines on the energy flux surface are those of AS(x), 
and hence indicate the direction plasma is moving. 



at distances r ~ 1.5 fT from the quark. Therefore, essentially all of the structure seen in Fig. [3] 
is simply that of hydrodynamics. The agreement between hydrodynamics and gauge/gravity du- 
ality in a strongly coupled SYM plasma bolsters the notion that one should be able to model 
accurately the wake produced by a high energy quark (or gluon) moving through a strongly 
coupled QCD plasma merely using hydrodynamics. 




Figure 4: Left — A log-log plot of the light quark stopping distance Ax as a function of total quark energy E for many 
falling strings with different initial conditions. Right — The instantaneous energy flux down a falling string normalized 
by its initial energy Eq. 

The penetration depth of a light quark (or gluon) was studied in Refs. |[T5"ll23l . In Ref. ||231 it 
was argued that the total distance a light quark can travel scales with the cube root of its energy in 
the high energy limit. Fig.|4j taken from Ref. [ 15], shows the penetration depth of nearly-onshell 
light quarks as a function of energy for many different sets of initial conditions and for energies 
E » V/LT^] As is evident from the figure, for a given energy the total distance a quark can travel 
is not unique. However, for a given energythere is a maximum distance a quark can travel. This 
bound, shown as the solid red line in Fig.Ul is given by Ax max = (0.526/T)(E/T V^) 1 ' 3 . 

From the E 1 ^ scaling of the penetration depth one might naively expect the light quark 
energy loss to scale like dE/dt ~ E 2 ^. However, as discussed in Ref. [ 15] this turns out not to be 



3 In the gravitational setup the different sets of initial conditions for a quark simply amount to different initial condi- 
tions for a falling string. 
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the case. To elucidate this point, also shown in Fig.[4]is the instantaneous energy flux down a dual 
falling string as a function of time. The energy flux is evaluated a spatial distance * 1/2T from 
the string endpoint. As discussed in the previous section, the flux of four momentum down the 
string is related to the rate the dual quark looses four momentum to the plasma. More precisely, 
the energy flux down the string — averaged over microscopic times, which are order l/T — 
yields the four vector f appearing in the effective source for hydrodynamics in Eq. d5l. As is 
evident from Fig.|4],the energy flux down the string does not decrease in a power-law fashion as 
a naive E 2 ^ scaling of dE/dt would suggest, but rather increases monotonically ! 

As stressed in Ref. [ 15], the precise form of the energy flux down the string is sensitive to the 
initial conditions used to create the string. However, for very long-lived excitations the late time 
behavior of the energy flux is universal. At late times the energy flux away from the string end- 
point grows like 1 / V f thei-m _ t where f t herm is the thermalization time {i.e. roughly the time where 
a falling string's endpoint approaches the event horizon) [15 |j^] This late-time behavior implies 
that after traveling substantial distances through the plasma, the thermalization of light quarks 
ends with a large transfer of energy to the plasma. This behavior is qualitatively similar to the 
energy loss rate of a fast charged particle moving through ordinary matter, where the energy loss 
rate has a pronounced peak (known as a "Bragg peak") near the stopping point. It is noteworthy 
to mention that this behavior is qualitatively different from that of heavy quarks. At both weak 
and strong coupling heavy quarks prefer to deposit most of their energy and momentum during 
the initial stages of their trajectory, when their velocity is largest Ifl2ll24ll . 
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